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We attempt to deal with the orbifold singularities in the moduli space of flat con¬ 
nections for supersymmetric gauge theories on the torus. At these singularities 
the energy gap in the transverse fluctuations vanishes and the resulting break¬ 
down of the adiabatic approximation is resolved by considering the full set of 
zero-momentum fields. These can not be defined globally, due to the problem of 
Gribov copies. For this reason we restrict the fields to the fundamental domain, 
containing no gauge copies, but requiring a boundary condition in field space. 
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This paper is dedicated to the memory of Michael Marinov. He 
and his family suffered the hardship and humiliation of a “refusenik” 
in the former Soviet Union. I admired him for his strong moral prin¬ 
ciples and persistence. I feel fortunate having known him, and having 
experienced his kindness. Michael has done pioneering work involv¬ 
ing Grassmann variables, supersymmetry, geometric quantization and 
quantum tunneling. I hope the following result would have been to 
his liking. My attempts to address this problem stem from the period 
I first met Michael at a workshop in Trieste, now just over 10 years 
ago. 


1 Introduction 

We revisit supersymmetric Yang-Mills theorieSj-pn the torus to study the vac¬ 
uum state in connection with the Witten indexS The torus geometry is crucial 
to preserve the supersymmetry. The index counts the number of quantum 
states (fermionic states with a negative sign). Due to supersymmetry, states 
at non-zero energy occur in fermionic and bosonic pairs, and do not contribute 
to the Witten index. The counting can therefore be reduced to the vacuum 
sector. The Hamiltonian is given hy H = ^{Q, Q^} with Q, the supersym¬ 
metry generators, and unbroken supersymmetry requires these to annihilate 
the vacuum state, hence giving a zero vacuum energy. Only at zero energy 
there can be an absence of full pairing, which would make the Witten index 
non-zero. A zero value of this index thus indicates supersymmetry may be 
spontaneously broken. 

In perturbation theory bosonic (gluon) loops are canceled by fermionic 
(gluino) loops, and applied to the problem of non-abelian gauge theories in a 
finite volume, it leads to the absence of an induced effective potential on the 
moduli space of flat connections (the so-called vacuum valley). The gluinos 
are represented as Weyl fermions in the adjoint representation of the gauge 
group, denoted by A“, with a a two-component spinor index. They are the 
superpartners of the gluons. 

A technical problem that has remained unresolved ever since Witten’s 
original work is associated to a breakdown of the adiabatic approximation in 
the reduction of the degrees of freedom to those of the classical vacuum, when 
using periodic boundary conditions. This classical vacuum is defined up to a 
gauge transformation by the set of zero-momentum abelian gauge helds. Its 
gauge invariant parametrization is in terms of the Wilson loops that wind 
around the three compact directions of the torus, which are compact variables. 
This describes the vacuum valley as an orbifold, T^/Z ’2 for SU(2). The orbifold 
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singularities arise where the flat connection is invariant under (part of) the 
Weyl group (the remnant gauge transformations that leave the set of zero- 
momentum abelian gauge fields invariant). For SU(2) their are eight orbifold 
singularities (related to A = 0 by anti-periodic gauge transformations). 

Without the contributions from the fermions, the wave function would be 
localized to these orbifold singularities. The singularity is resolved by includ¬ 
ing all the zero-momentum gauge fieldsoa A reduction near A = 0 in terms 
of the abelian modes is impossible, since here the energy gap between the 
fluctuations in the abelian and the non-abelian held directions vanishes. This 
is the source of singular non-adiabatic behavior, and remains so for the su¬ 
persymmetric caseQ Resolving the orbifold singularity in the supersKmmetric 
case leads one to consider the reduction to zero (spatial) dimensionsBQ In the 
context of the supermembraneQ it was found that the spectrum is continuous, 
down to zercfaenergy. One can construct trial wave functions with arbitrarily 
small energyP by moving its support far from A = 0. In the case of the su¬ 
permembrane the vacuum valley is non-compact. For gauge theories on the 
torus, the compactness of the vacuum valley arises due to identifications under 
periodic gauge transformations with non-zero momentum. Such a gauge trans¬ 
formation does not preserve the momentum for the non-abelian modes. What 
is zero-momentum near A = 0, is non-zero momentum near a gauge copy of 
A = 0. We therefore have to match, somehow, the behavior near each of the 
orbifold singularities to the behavior far away, where a reduction to the abelian 
zero-momentum modes is dynamically justified, and where the Hamiltonian is 
just the standard Laplacian on the torus. 

1.1 Status quo 

Let us first summarize the current state of affairs. Assuming the reduction to 
the vacuum valley is justified, we note that the zero-momentum gluinos asso¬ 
ciated with the abelian generators, each witlutwo helicities, carry no energy, 
which is the source of the vacuum degeneracylll These gluinos have to be com¬ 
bined in Weyl invariant combinations, respecting Fermi-Dirac statistics. There 
are r (with r the rank of the gauge group) independent invariants, made from 

U = (1) 

and its powers. So one has C/"|0), n = 0,1, • • • ,r, as bosonic vacuum states, 
and no invariant fermionic vacuum states. This led Witten to an index equal to 
the rank of gauge group pluSpOne, r -I- 1. To circumvent the problems near the 
orbifold singularities, Wittenlll considered the alternative of twisted boundary 
conditionsu For SU(IV) the same result for the index follows. However, other 
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groups do in general not admit the type of twisted boundary conditions that 
completely remove the continuous vacuum degeneracy (with its associated orb- 
ifold singularities). In particular these twisted boundary conditions could not 
be used to resolve a discrepancy with the infinite volume result based o n th e 
determination of the gluino condensate through instanton contributionsliJ’tiiEj 
relying on the fact that the index should not change with the volume, or for 
that matter any other smooth deformation. 

The gluino condensate calculation has been justified by first adding matter 
fields, which introduce an external mass scale so as to control the weak cou¬ 
pling expansion. One then relies on the index being constant under a smoo^ 
deformation (through holomorphy), that decouples the extra matter sectortH 
In the direct (strong coupling) approach, since the instantons have more than 
two gluino zero modes, it seems the condensate {XX) vanishes. Instead the 
appropriate power of the gluino condensate, ((AA)^), is considered where h 
is the so-called dual Coxeter number of the gauge group, h = N for SU(fV), 
which counts the number {2h) of gluino zero modes. Invoking the cluster de¬ 
composition property, it is this power in h that gives the number of vacuum 
states, 

(AA)=e2-“//^(|((AA)^)|)'/^ n = (2) 

These arguments seem reasonable, but are not rigorous^lli and suffer from a 
discrepancy in the prefactor (-^5/4 for SU(2)). More recently, u^lms been 
made of the constituent nature of periodic instantons (or calorons)li3’t3 in the 
context of a Kaluza-Klein reduction with periodic gluinos (as opposed to a 
high temperature reduction with anti-periodic gluinos, which would break the 
supersymmetry). The constituent monopoles, with Aq playing the role of a 
Higgs field, have exactly twi zero-modes and saturate the condensate, (AA), 
giving the correct prefactorE3 In this case it is the compactification scale that 
justifies the semiclassical approximation. 

The mismatch in the Witten index between small and infinite volumes 
occurs for SO{N > 6) and the exceptional groups. There has, however, been 
a recent revision in counting the number of vacuum states ixL a finite volume. 
In a study of D-brane orientifolds in string theory, WittenU constructed for 
SO(7) an extra disconnected component on the moduli space of flat gauge con¬ 
nections, which can be embedded easily in SO(A^ > 7). For SO(7) and SO(8) 
this gives an isolated component of the moduli space, contributing only one 
extra vacuum state. For SO{N > 8) the extra component in the moduli space 
behaves like the trivial component for SO(W7). Adding r-|-1 coming from the 
SO{N) and SO(A^— 7) moduli space components gives the dual Coxeter num¬ 
ber of SO{N), thereby yielding the same number of vacuum states as obtained 
in the infinite volume. 
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Witten’s construction based on orientifolds does not work for the excep¬ 
tional groups. This natr^lly led to a derivation of the extra vacuum states in 
a field theoretic contextli3 trivially extended to the exceptional group G 2 , as 
a subgroup of SO(7). It was subsequently solved for other exceptional groups 
with periodic baundary conditionsEjCu and for any group with twisted bound¬ 
ary conditionsEj Twisted boundary conditions usually do not remove all the 
vacuum degeneracies, but it is important that the number of vacuum states 
is independent of the twist for all gauge groups that have a non-trivial cen¬ 
ter. The origiiuof the extra moduli space components is actually not too hard 
to understandllj Large gauge groups can have subgroups that are products 
of unitary groups, which each would allow for twisted boundary conditions. 
By choosing “twists” from all subgroups to cancel (or give the desired total 
“twist”), one obtains flat connections that can not be deformed to the Cartan 
subalgebra (which supports the trivial component of flat connections). 

Although these new results for counting the number of vacuum states in a 
finite volume remove the urgency of addressing the problem with the adiabatic 
approximation, it does r|emain a sore point in the finite volume analysis, as 
stressed again by WittenB3 On the one hand the wave function on the vacuum 
valley has to be constant, on the other hand it seems to want to vanish near 
the or bifold singularities. We will argue that deviations from a constant will be 
confined to a distance from each orbifold singularity that is 0{g^/^{L)) times 
the distance between the orbifold singularities, where the dependence on L is 
due to the asymptotically free non-trivial running of the coupling, appropriate 
for = 1 supersymmetric gauge theories. At the same time, however, in this 
small region the energy needs to remain zero and here lies the burden of the 
proof. 

1.2 The adiabatic approximation 

We are interested in constructing the vacuum wave function in sufficiently small 
volumes. Our convention is to choose the dependence on the bare coupling 
constant such that it appears as an overall factor I/pq, 

= + (3) 

The reduction to the zero-momentum degrees of freedom, as in the bosonic 
case, will replace the bare coupling constant by a running and asymptot¬ 
ically free coupling constant g{L). The zero-momentum gauge fields are 
parametrized as Ai = ic'fTa/ (2L), with t the Pauli matrices. The vacuum 
valley is parametrized by the abelian degrees of freedom. These are defined by 
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Ti, with TiTj = for each i and j. Alternatively, we may parametrize 

the vacuum valley by = Q = cf, choosing the maximal abelian subgroup 
to be generated by the third Pauli matrix T 3 . The effect of the periodic gauge 
transformations, 

gn(x) = exp(-27rm • xr^jL), (4) 

is to shift C over 47 rn, making the vacuum valley into a torus, see Fig. 



Figure 1: A two dimensional slice of the vacuum valley along the (Ci,C 2 ) plane. The grey 
square is the fundamental domain A. The dots are gauge copies of the origin (which turn 
out to lie on the Gribov horizon fl, indicated by the fat square). 

Also in the supersymmetric case, the Hamiltonian is invariant under anti- 
periodic gauge transformations (gauge transformations periodic up to an el¬ 
ement of the center of the gauge group). When n has at least one of its 
components half integer it is homotopically non-trivial. Since it is a symmetry 
of the full Hamiltonian, this is one way to see that the vacuum valley has 8 
orbifold singularities, related to A = 0 by shifts over half-periods (27r in each 
of the three directions). Alternatively, the action of the Weyl group on the 
vacuum valley is given by the reflection C ^ —C, and the orbifold singular¬ 
ities are at the eight fixed points of this symmetry (combined with the shift 
symmetries). The cell C G [— 7 r, 7 r]^ can be used as a fundamental domain A, 
see Fig. |l|. Any point on the vacuum valley can be reached by applying suitable 
gauge transformations. Opposite sides on it’s boundary are identified under 
these homotopically non-trivial gauge transformations. The representations of 
their hpmotopy define the electric flux quantum numbers as introduced by’t 
Hooftcl We will here only consider the sector with zero electric flux, i.e. the 
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trivial representation, where wave functions at opposite sides are equal. 

We need to reconsider the construction of the effective Hamiltonian, since 
the gluino loops tend to cancel the gauge loops. For a background with zero 
field strength, the effective potential is easily seen to vanish to all orders in 
perturbation theoryH But to resolve the orbifold singularity near A = 0, we 
do not wish to integrate out the non-abelian zero-momentum modes. These 
modes can have non-zero field strength, and the quantum corrections for these 
are expected not all to cancel. Otherwise the /3-function would vanish, which 
we know not to be the case for = 1 supersymmetric gauge theories. If 
necessary, field redefinitions should restore the supersymmetry. 

In the background field gauge the one-loop effective potential reads 


-1 


1 
4 


g 2 — o: 2 {d) J (F'jy) -l-a 3 (F“) (cj.) -I-osdet (c)-I-■ 


Vi{c)=L 

l u V .(/n / " " J 

J 

where = —EabdC^jcf.- The coefficients (labeled as in the bosonic casd23) 
in dimensional regularization are, up to terms vanishing at d = 3, given by 
(fc G (27rS)3) 

/ j\ {d 3)(d -)- 6 ) \ 1 1 \ 1 15 \ 1 

“ 12 d ~ 

k^O ' ' fe/O ' ' fc#0 ' ' 


.( 6 ) 

The result for dimensional reduction (vector and spinor indices strictly in 4 
space-time dimensions) is simply obtained here (and below) by putting d = 3, 
oii 02 ( 3 ) = -f EfeVo Fl”' = 4 ^^ + ^ (it + 0.409052802 • • •) +0(d-3). 
The effective potential vanishes along the vacuum valley, as expected. 

In minimal subtraction one defines 


- 47rHd-3)^ g^iLY g^L) - 4 ^^^^ ms) 


( 7 ) 


such that L‘^~^jgl + a 2 (d) = lfg‘^{L) + with 02 the finite part of Qf 2 (d), 
and (/(T) the running coupling appropriate for the supersymmetric theory. The 
“electric” part of the effective Lagrangian is also not difficult to compute in 
the background field calculation, 


ML 


d-3 


9o 


+ ai{d) ] {{DodY{t)Y , ai(d) = - 


3-l-5d^ 1 

4d 

k^O ' ' 


( 8 ) 


With {DodY = cf + EadbC^ci, we keep eg ^ 0 to preserve supersymmetry, 


^7l“ = -(e- 7 MA“-A“ 7 ^£), dA“=7'^%V, dA“ =-£ 7 ^"^“ , (9) 
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which allows for a consistent truncation to the zero-momentum sector. Both 
in dimensional regularization and dimensional reduction the finite parts, ai^ 2 i 
are equal and will be absorbed in the running coupling. Finally, for the gluino 
part of the effective Lagrangian we find 


iZ 


(^F7oA“(t) -I- K£adbcf{t)X^{t)^ 


Z = l 


9q 


E 

fe #0 


|fc|3 


K, = 1. 


( 10 ) 

In the background held gauge only the gluino Held requires a rescaling to 
restore the supersymmetry (with our use of component Helds, implying the 
Wess-Zumino gauge, we do not maintain explicit supersymmetry, but the final 
results should of course not depend on this). j—. 

We also used the Lorentz gauge, which shifts the coefficients byEa 
S{ai,a 2 } = {l, 2 }X)fe 5 .io d{a3, Sa^} = ■ In ad¬ 

dition 5 Z 2 = 0, and 5k = ^ai. Infinities are absorbed by a simple rescaling of 
the gluon field, as guaranteed by gauge invariance. The changes in as and as 
can be absorbed by a finite but non-luifiar field redefinitionCJ All we want to 
stress here is, like in the abelian case^ that supersymmetry does not imply 
that the corrections (like as and as) need to vanish, but a self-consistent cal¬ 
culation requires further terms in the effective action, involving a non-static 
backgrpund field (for the bosonic case most reliably done in a Hamiltonian 
setting23). The importance of these higher order terms will become clear later. 
To lowest order, one finds the Hamiltonian truncated to the zero-momentum 
sector, with the bare coupling replaced by the renormalized coupling. 

The energy gap in the fluctuations transverse to the vacuum valley is easily 
read off from the Lagrangian. Close to the origin it is given by 2\C\/L. Inte¬ 
grating out transverse degrees of freedom is only reliable if the energy of the 
low-lying states is smaller than this gap. This energy behaves as g‘^^^{L)jL 
(as shown in the next section, rescaling c with makes the Hamilto¬ 

nian proportional to g^^^{L)lL). The gap in the transverse fluctuations thus 
becomes of this order where \C\ is of order 

Consider a sphere of radius g^/^[L) around each orbifold singularities, be¬ 
yond which the adiabatic approximation is accurate. There, the wave function 
can be reduced to the vacuum valley and, assuming indeed it has zero energy, 
it will be constant. (Note that the chosen parametrization of the vacuum val¬ 
ley is independent of g{L) and L.) As long as g{L) is small, we may assume 
the wave function in the neighborhood of the orbifold singularity to respect 
spherical symmetry. In the bosonic case, once the wave function spreads out 
over the vacuum valley, any spherical symmetry is quickly lost. In the super- 
symmetric case, however, the reduced wave function on the vacuum valley will 
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become constant up to exponential corrections at distances much greater than 
from the orbifold singularities, which are separated over a distance 2 tt. 
(Rescaling c with the boundary of the sphere introduced above is at 

g~^/^{L), with the boundary of the fundamental domain and the other orbifold 
singularities an other factor g~^/^{L) removed from that.) Instead of insisting 
that the groundstate wave function is normalizable, we should rather insist on 
its projection to the vacuum valley to become constant. As Jang as the wave 
function is bounded everywhere, the^issues of normalizabilit'^j in the context 
of the supermembrane or M-theoryCa applications, is therefore of no relevance 
here. But the challenge of explicitly constructing the ground state, of course, 
remains the same (a subtle difference will be pointed out later). 

As we have argued, when the volume is small enough, we may assume 
the groundstate wave function of the effective Hamiltonian to be spherically 
symmetric. The boundary conditions at the boundary of the fundamental 
domain are replaced by requiring this wave function to approach a constant, 
after projecting to the vacuum valley. This projection is well defined, because 
at large separations from the origin, the wave function becomes exponentially 
localized transverse to the vacuum valley. Spherical symmetry near the orbifold 
singularities will dramatically simplify the analysis. 

In the next section we will set up the zero-momentum supersymmetric 
Hamiltonian, and its rediiction to the gauge invariant, spherically symmetric 
sector. This is not new^GS but we will be able to push it to the point where 
we can explicitly construct a complete basis of states that respect these sym¬ 
metries. The Hamiltonian can be split into a “radial” and “angular” part, and 
our basis explicitly diagonalizes the angular part (“spherical harmonics”) in 
terms of invariant polynomials. This may be useful in a more general context. 


2 The Hamiltonian 

The conventions for the superalgebra we follow are those of Wess and BaggerS 
Jie supersymmetric Hamiltonian we follow to a large 
i,po| yyg start from the supercharge operators, 


For the formulationn£ 
extent earlier workE3E 


Qa — 







( 11 ) 


with = cf/{g{L)L) and (and the unit) as 2 x 2 matrices. Re¬ 

stricting to the zero-momentum modes, both the Weyl spinors Af and Af 
are constant. Lowering indices is done with e^/? = = —iT 2 , 6ab and 

= diag(l, 1,1, —1) (or Sij) for respectively the spinor, group and space- 
time (or space) indices, and raising of indices is done with the inverse of these 
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matrices. Repeated indices are assumed to be summed over, but to keep no¬ 
tations transparent we will not always balance the positions of the gauge and 
space indices. For zero-momentum gauge fields 

= -^.gSrjkSabcV^V^. ( 12 ) 

In the Hamiltonian formulation the anti-commutation relations 

{^aa^^b/3} ^ {A““,A^^}=0, {A““,P^}=0, (13) 

with (To the unit 2x2 matrix (one has ((t^)““ = e“^e“^((T^)^^), give 


{Qc, Qa} = 2((To)aaH - 2((T*)„^y“^„ (14) 


where 

Ga = igsabc - A^doA^ j (15) 

is the generator of infinitesimal gauge transformations, and Ti. is the Hamilto¬ 
nian density 


H = -i 




dV‘^dV‘^ 


+ - igeabcX‘^a^X\^. 


(16) 


Splitting the Hamiltonian, f d?xH = g‘^/^{L)H/L, in its bosonic and fermionic 
pieces, H = Hb + Hf, we find with c“ = g'^^^{L)c°; 


+i{Bty, Hf = -ieabdX^a^X^cf, (17) 

where 

K = -le^^'^eabdc’^ct. (18) 

As discussed in the previous section, the orbifold singularities, other than 
at c = 0, lie at a distance 27rg“^/^(L) in these new variables c (measured along 
the vacuum valley where B vanishes). We want to solve for the groundstate 
wave function such that for |c| ;g> 1 it becomes a constant, after projecting 
to the vacuum valley (we will come back to this projection later). As this 
boundary condition is compatible with spherical symmetry, i.e. it goes to the 
same constant for all directions on the vacuum valley, we will restrict ourselves 
to wave functions ^'(c) that are spherically symmetric and gauge invariant. 
We stress this is an accidental spherical symmetry, that holds in sufficiently 
small volumes. 



The Witten index 11 


Building the Fock space of invariant states, we first separate in the fermion 
number .|_Sates with odd fermion number do not respect the symmetry and we 
can onlycZl have F" = 0, 2,4 and 6 (there are six independent Weyl components). 
Particle-hole symmetry relates F = 0 to F = 6 and F = 2 to F = 4. Since 
Hf\F = 0) = Hf\F = 6) = 0 (the diagonal entries of it* vanish), this case 
reduces to the bosonic Hamiltonian which is known not to have a zero vacuum 
energyB So we can restrict our attention to the |F = 2) states. The index will 
be twice the number of zero-energy states in this sector (due to the particle- 
hole symmetry, see below). It should be noted that only the relative fermion 
number is well-defined, because integrating out certain modes involves filling 
negative energy (one-particle) gluino states. 

2.1 Invariant two-gluino states 

Instead of making irreducible decompositions of the variables^ we chose to 
write down the most general F = 2 states, and show how the Hamiltonian 
acts on these. We can combine two-spinors symmetric or antisymmetric in the 
gauge index (and thus respectively antisymmetric and symmetric in the spinor 
index) 


|V) = = -2iV;eahcA“(a^°)“P^|0), 

\S)=SabJ‘^’’ = -Sab~Xl~Xy^\0), (19) 


where — a^a^), such that as a 2 x 2 matrix with the 

first index up and the second down (bringing also the first index down leads to 
a symmetric 2x2 matrix). Here 12“ and 5“^ = are arbitrary (and assumed 
to depend on c“). The action of the Hamiltonian on these states preserves this 
structure 

Hf{\S) + |V)) = |V) + |5) (20) 

with 

V.“ = e,,keabdc]vt + - c\S^\ 

gab ^ 25‘^^cfvt - - c,^V“. (21) 

Covariance allows us to write the most general form 


V, 


5 


ab 


hi(f,u,v)cyf - /i 2 (f,zi,u)F“/f^ -t h 2 ,{f,u,v)c]c\cl/r^, 
hiir, u, u)(5“^ -f /i 5 (f, u, + he,{r, u, u)c“Cj c^Cfc/f"^, (22) 
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in terms of the invariants = (c“)^, u = and v = f ^detc. It is 

useful to introduce the polar decomposition!^ c“ = J2j with R,T G 

5'0(3), since the spherical symmetry and gauge invariance allows us to put 
c“ = diag(xi, a: 2 , a; 3 ). Note that this does not completely fix the freedom under 
gauge and spatial rotations. The remnant symmetry involves permutations of 
the Xi and simultaneously flipping two of its signs. In terms of x, 
u = ^i>j ^ rij > properly invariant under the remnant 

symmetry. Fixing this remnant symmetry, for example by |xi| < X 2 < X 3 , 
allows one to solve the Xi from {f,u,v). That V and S can be expanded each 
in terms of three invariant functions (the hm) is most easily established in this 
diagonal representation 

Vj“ = Sj (hi{f,u,v)f~^x'^ + h 2 {f,u,v)f~'^detc/xj + h^{f,u,v)f~^xf) , 

gab ^ ^ab _|_ h^(r, u, v)r~‘^x\ + he{f, u, v)f~'^xl) , (23) 


(note that Bj = —S^xJ^detc] no summations over repeated indices). Any 
higher order term can be reduced to this form. This is best illustrated by 
examples: one brings and to the respective form of V“ and 5“** 

in Eq. (p^), by using the identities Xj = f'^x^ — r^uxj + jxj and x\ = 


-2 4 
r xl 


-4 2 

— r ux„ 


fV. 


We have established that any invariant state |d>) can be decomposed as 


6 

1^) = X! h„,{f,u,v)\e„,{u,v)), (24) 

m—1 

where the |em) are implicitly defined by Eqs. ( |l^ , p^ . Since Hf does not 
contain any derivatives with respect to c, we can diagonalize Hf pointwize. 
First determine the matrix of Hf with respect to the basis |em) 




From Eq. (^) one directly reads off that 


HJ^^ = f 


( 0 1 
2 0 
0 -1 


— V 
1 
0 


2 1 l-u\ 

0 0 —u 

0 -1 -1 


2 4u 2 - 4 m 0 0 0 

-2 0 0 0 0 0 

VOO -2 00 0/ 


(25) 


(26) 
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This matrix is not symmetric due to the fact that the \em) are in general not 
orthogonal. The non-diagonal norm matrix = (em|e„) is however block 
diagonal, since (5|V) = 0 for any choice of V and S. Introducing the notation 


X = l-2u,Y = l-3u + 3v^, 

z = 

1 -t- 2 'u 

— Au + Av 

^ one 

finds 

/ 8 

2 Av 

8 X 

0 

0 

0 \ 


2 Av 

8 u 

8 v 

0 

0 

0 


j\^mn _ 

8 X 

8 v 

8 Y 

0 

0 

0 



0 

0 

0 

12 

4 

AX 



0 

0 

0 

4 

AX 

AY 



\ 0 

0 

0 

AX 

AY 

AZ ) 



The matrix N can be used to make an orthonormal basis. It transforms Hf to 
Ni ■HfN~ 2 . This is seen to be symmetric, e.g. by establishing the symmetry 
ot N -Hf. 



-1 - 0.5 0 0.5 1 

Figure 2: The six eigenvalues Efir, of Hfff, plotted versus 3u\/3. 

We denote hy Ef the energies of the two-gluino state in a given background. 
With Ef linear in f and, as we will see, independent of u, only the v dependence 
is non-trivial. The eigenvalues ot Hf, for the six invariant two-gluino states, 
are determined by 

det{Hf - Ef) = {Ej-Ac^Ef - 16detc){Ej - c^Ef + 2detc) = 0, (28) 

In terms of the three roots pa, E^—c^Ef+2det c = {E f—p,i){Ef—fi 2 ){E f—, 
the remaining three roots are given by pi = —2pi. Note that ^^Pi = 0, 
Ti = 2c2 and 0, Pi = — 2detc. Explicit expressions for the two-gluino 
energies are (see Fig. ||) 

Pk = —2pk = —4fcos(27rfc/3-I-arccos[—3u-\/3]/3)/'\/3, 


(29) 
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with arccos ranging from 0 to tt. Note that Sv'/S € [—1,1] since Yl^i 
extremal (under the constraint = 1) when all |a;i| are equal (to l/-\/3)- 

The lowest eigenvalue, = /is, takes on the value —2f along the vacuum 

valley, since u = 0 implies v = 0. The eigenvectors associated to the eigenvalue 
jlk are also easily constructed, 

= (i/ife/f, 1,0, l + 4wf/(ifc,-l,0)/|4(l + 6uf//lfe)|, (30) 

normalized in a fixed background according to rt=i ^ = ^kk'- 

We will not write down the more complicated analytic form for the eigen func¬ 
tions belonging to the remaining three eigenvalues. Only will be referred 
to later when discussing projection to the vacuum valley. 

At this point it is perhaps interesting to observe that along the vacuum 
valley (two out of the three Xj vanish), the particle-hole dual of |'k) = hm\em) 
can be expressed as |^) = ^A[(,c“c(’A“|4'), obtained from |4>) by creating a 
pair of_abelian zero-momentum gluinos, in accordance with the analysis of 
WittenEJ To prove this, note that with c“ abelian, defining V“ = —cf and 
Sab = Sab — c“Ci, we Can write Idt) = |V) -I- |5). Its particle-hole dual is given 
by 1^) = |V) -I- |5), with 

|V) = = -2zV;ea6cA“(a^°)^A'^|0), 

|5)=5a6cJ“'> = -5,feA“A^e“^|0), |0)=nA“|0), (31) 

acx 

also defining the particle-hole symmetry in the general case. 

It should be noted that our change of variables from c to Xi will suffer from 
coordinate singularities, also affecting the linear independence of the jem). This 
is most obvious from the fact that det N = 32768with J (up to a constant) 
the Jacobian for the change of variables, d^c = f tt'* (with jxi | < a ;2 < 2 : 3 ), 


J = '^\xl — x^j \ = — Au) — z)^(4 — 18m -I- 27z)^). 


(32) 


i>j 


Indeed, the normalization of in Eq. (30) is singular when fik + Gvf vanishes, 
which is easily identified with a vanishing Jacobian. For this observe that 
ic^Ef J- 18 det c = + 6vf) equals the difference of Ej — A<?Ef — 16 det c 

and Ej — (?Ef + 2detc. Since each of these vanish for either Ef = fii or 
Ef = p-i, we conclude that fik + 6nf vanishes if and only if fik coincides with 
one of the other roots Hi. From Fig. || we read off this happens for v = ±-\/3/9, 
which is the extremal value of v where all three |xi| are equal (constraining u 
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to be 
exists, e.g. 


1/3), and therefore J = 0. As long as J ^ 0, the limit w ^ 0 


lim 

11^0 




i 0 i -1 0\ 
100 - 10 , 
i 0 1 -1 oy 


(33) 


independent of u, but we have to remember that when J = 0 (in particular 
M = 0) the I Cm) are no longer independent. We come back to this when 
projecting to the vacuum valley. 

Before constructing a full invariant basis, to gain more confidence that the 
I Cm) correctly describe the invariant two-gluino states, it is useful to note that 
the spectrum oi Hf can also be obtained from its (non-invariant) one-particle 
statest3 For this we write 


Hf = (34) 

The one-particle fermion energies are given by the eigenvalues of M. Note 
that charge conjugation symmetry, M* = T 2 MT 2 , implies each eigenvalue is 
two-fold degenerate. By direct computation one verifies thatt3 

-c^M -2Aeic = Q, (35) 

such that any eigenvalue /i of M satisfies the equation /j,^ — — 2 det c = 

+ fL 2 ){T + Ts) = 0) with fii as defined before. Hence the one-particle 
fermion energies are given by —fH: and each occurs with a two-fold (spin) 
degeneracy. We can make invariant two-particle states only by combining two 
opposite spin states. There are six such invariant states, three with the same 
one-particle energies, giving Ef = —2fLi = p,i, and three with different one- 
particle energies, giving Ef = —fii — fLj = Hk with i ^ j ^ k. This agrees 
with our earlier results. Note that there are also three triplet two-gluino states 
with two-particle energies Ef = fik, combining two one-particle states with 
different energies but equal spin, which are not invariant and thus left out 
from our considerations. See the Appendix for yet another method. 


2.2 Supersymmetric spherical harmonics 


We have removed the angular degrees of freedom on which the invariant wave 
functions do not depend: those associated with the gauge and space rotations. 
This leaves a six component F = 2 wave function., depending on either (f, m, v) 
or equivalently {xi,X 2 ,X 3 ). It is well knownrnrn that the kinetic term of the 
Hamiltonian reduces to 


52 

Hdcfr 


d 


d 


-iJ 


dxj 


dxi 


,_8_5 

df^ dr 




(36) 




16 P. van Baal 


with A{u,v) the Laplacian for S'®, acting on invariant functions, given bjt3 

q‘2, q‘2, q‘2 , 

A{u,v) = 4i3v‘^ + u-4u^)j^ + 8il-3u)v-^ + {u-9v^)j^ 

+4(2-1I„)|--3o4. (37) 

Note that for the bosonic theory {F = 0) the variables (f, u, v'^) were used (for 
the positive parity states), but for the F = 2 sector functions even and odd in 
V mix, requiring us to use the variables {f,u^v). The Hamiltonian still allows 
for a radial decomposition (note that df\em) = 0) 

or or 

where see Eq. (||), and 

6 6 

\en{u,v))H'^hm{r,u,v) =-^A{u,v)'^hm{f,u,v)\em.{u,v)). (39) 

m,n=l m—1 


It is straightforward, but quite tedious, to explicitly calculate 

= (40) 

For the first term A only acts on the coefficient functions, for the next we 
collect the terms where one derivative acts on the coefficient functions and one 
on the two-gluino basis vectors \em), whereas the last term has all derivatives 
acting on these basis vectors. This splits in the two sectors associated to |V) 
and |5), since derivatives cannot lead to mixing of these two-gluino states, 
specified by X and J in Eq. ( |T^ . We find the following results 

/ (2 — ^u)du — 3vdv dy + 2vdu 3vdy + 6udu \ 

Ay = 2 dy (2 - 8u)du - 6vdy -6vdu , (41) 

y —2du —dy —12udu — 9vdyj 

/O 2vdy 

A ;5 = 2(0 (4 — 8u)du — 6vdy 

Vo -Adu 

and 

/-4 0 7 \ 

A° = 2 0 -9 0 , A;^ 

Vo 0-15/ 


-8v‘^du 
4vdy + 8udu 

-IQudu — 12vdy 


(42) 



( 43 ) 
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Invariant wave functions, 
ized using 

(^1^') = / 


= Yfm.=i^rn{f,u,v)\eni{u,v)}, are normal- 
d^c ^ (44) 

m,n=l 


with iV"*” the norm matrix as defined in Eq. ( p^ . The matrix elements of the 
Hamiltonian with respect to such a basis are thus given by 


(4'|i7|4'') 



(45) 


with the matrix operator as defined in Eq. (p^. 

Like for the bosonic case (with H/\ = — iA and 77/ = 0) we first construct 
an invariant set of “angular” wave functions as polynomials in u and v. These 
can be chosen to be eigenfunctions (3^f^^) of 77 a, proportional to coefficient 
functions defined by (for F = 0 : z = 0; s € {2, 5} and F = 2 : z G {1,2,3}; 
s €{0,1,2, 3}) 

/U^\ /O 0 0 1 0 0\ 

{72 = 0 0 0 0 1 0 

yt/oV yo 0 0 0 0 1 / 

/u^\ /o 1 0 0 0 o\ 

7/| = u 0 0 0 0 0 

\uiJ yo 0 V 0 0 0 / 

judiciously chosen such that q = 2x + 3y + z allows us to order states, and 

= i(2g + s - 2)(29 + s + 5)u^v^yU^ + (47) 


/t/A /l 0 0 0 0 0\ 

T/f = 0 0 1 0 0 0 , 

\ufJ yo u 0 0 0 0 / 

Ui = vU§, U° = vU° = v, (46) 


with SR^y^ a linear combination of these monomials, but of lower order in q 
(states with different s do not mix under 77a, but they do mix under Hf). 
The appearance of the combination 2x + 3y is related to the fact that u^v'^y = 
^-2{2 x+d.y){ jIqI cj'^y _ The eigenvalues of 77 a can be written as L{2L + 7), 
with the “angular momentum” L = (7 + s/2—l = 2a; + 3?/ + z + s/2—1 taking 
half integer values. Eq. (|^ allows us to solve [77a - L{2L + 7)\y^y^ = 0 for 
j^xyz order g < 2a; + 3^ + z), with 

y^y^ = u^v^yu^ + R^.yy ( 48 ) 


To order these states at degenerate values of L, we begin at 7 = 0 with 
yi = 3^0°^ ■ Raising L with i we start with the lowest value of q and construct 
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Table 1: All orthonormal spherical harmonics for L < 4, with L = 2a: + 3y + s/2 — 1, 
such that = L{2L + For higher L use the available programcSl 


L 

yn =y^^^={hl,h2,h3,h4,h5,he) 

0 

j'l = = ( 0 , 0 , 0 , 1 , 0 , 0 ) v'35/2/8^^ 

1 

2 

3>2 =3>ri = ( 1 , 0 , 0 , 0 , 0 ,0)yTM/167r2 

1 

3>3 =3>2°°^ = (0,1, 0,0,0,0)yTT^/167r2 

1 

^^4 =3>o°°" = ( 0 , 0 , 0 ,-i,l,0)3v^/167r2 

3 

2 

3)5 = 3)3001 = (0,0, 0, V, 0, 0) VT5M5/167r2 

3 

2 

J'e =3>r = (-^, 0 , 1 , 0 , 0 , 0)11 v47^/32^2 

2 

3)7 =3)o° 2 = 0,0,0,0)39v^/327r2 

2 

^8 = yr = ( 0 , 0 , 0 , -4i, l)429y77M/167r2 

2 

J'g = = ( 0 , 0 , 0 , + u, ||)3v'6149/2/16^^ 

5 

2 

3^10 = 3^3^"^ = (0,0,0, -f, V, 0)3v'51051/2/16^^ 

5 

2 

3^11 = 3>r^ = (-u, 0,0,0)39^/17^/64^2 

5 

2 

3^12 = yr = (-3 + u, ^, 0 , 0 ,0)33y2W7/167r2 

3 

3^13 = 3>2°°" = ^^, 1 ^, 0 , 0 , 0)51 v/4M^/327r2 

3 

3^14 =3>r = (-f|> -i + 11 :1': 0,0,0)39v'17765/7/647r2 

3 

3>15 = 3>o°“ = (0,0, 0, gig - ^ + u2, 0, 0)663y20977/327r2 

3 

3^16 = 3^1''" = ( 0 , 0 , 0 , ^ - ff, + u, ^)5 V2717/7/327r2 

7 

2 

3>i7 = 3>3°°" = ( 0 , 0 , 0 , if,-^, u)969v/23T7T0/32^2 

7 

2 

3^18 = 3 ^ 3 '''' = (0,0, O, -if + uv,-^, if )39\/53295/2/327r2 

7 

2 

3^19 = 3>1°“ = (^ - tI + -ff: - 353 : 0: 0: 0)969^42^/32^2 

7 

2 

3>20 = yr = (^ - It + ¥: + U, 0,0, 0)19^51051/2/64^2 


the independent, but in general non-orthogonal states Using that y and 

2 are uniquely fixed by g — 2 x, we start with a; = 0 and every time we increase 
X by one, we modify by projecting it on the orthogonal complement 

of the previous states (the Gramm-Schmidt procedure). When completed we 
increase q until all states with a given value of L are constructed, after which we 
increase L. Dividing by the norm we thus obtain a complete orthonormal set 
of “spherical harmonics”, yn{u,v) = {u,v\n), with n labelling the n-th state 
thus constructed. The first few spherical harmonics yn{u,v) are collected in 
Table |[ Note that to evaluate inner products we need to compute^e integrals 
For this we recall the recursive definition^ 

Y _ 4i(l + z + 4:j)Xi-ij + 12i{i — l)Xi-2,j+i + 2 j( 2 j — l)W+ij-i , 

(4i + 6j)(4i + 6j + 7) 
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with = 327r'^/105. The Mathemadcaa code for generating the is 
available through the World Wide WebE3 

We denote by the value ofL = 2a: + 3?/ + z + s/2—1 implicitly defined 
by yn = It will be convenient to also introduce £ = 2L + 3. With the 

spherical harmonics we can now construct a reduced Hamiltonian, 

{n'\H\n) = 77” ”(r) = K{r; 7v„)i5" " + ^f^{n\u\n) + r(n'|77/|n), 

K{f-L)=-if-^dff^df + f-^L{2L + 7) = r-^k{f,£)f^, 

k{r] i) = —^r~'^dfkdf + + 1). (50) 

We illustrate its sparse nature by showing in Fig. || the entries where either 
(n'|M|n) ^ 0 or (n'|77/|n) ^ 0 as black squares. The different bands can be 
traced to come from the selection rules |i5L| = 0,1,2 for the matrix elements 
of u, and \SL\ = I for the matrix elements of 77/. The number of codiagonals 
is bound by — 3 + 4y/n. 



Figure 3: Band structure of the reduced Hamiltonian, for the first 100 states. 


To write down the matrix of the full Hamiltonian with respect to an in¬ 
variant basis, we introduce radial wave functions (t>p{r) = {r\p,£). The radial 
quantum number p is associated to a momentum if we choose 


4>pik = CJ f) 


i\2 


K = m) 


(51) 
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where ji{z) is the spherical Bessel function which satisfies the equation 
K{f;£)ji{f) = or {f). Normalization 

with a factor Cp is such that / df r^(j)p{f)*(j)g{f) = 6pq. In terms of the nor¬ 
malized basis of invariant states (c|p, n) = (m, u|n)(f|p, €„), the matrix of the 
Hamiltonian is given by = (p',n\H\p,n), or 

^n'n =Ep^6nn'5PP + {p\ in'\r‘^\p, in) {n'\‘^\n) + (p'Jn'lflp, In) {u'\Hf\n) . (52) 

To complete the construction of the basis, we need to address the question of 
boundary conditions (fixing the momenta), which somehow have to incorporate 
that the groundstate wave function becomes a constant, after projecting to the 
vacuum valley. 

3 Vacuum valley and boundary condition 

The vacuum valley is characterized by those configurations for which u = 0 
(this implies v = 0), and f measures the (three-dimensional) distance to the 
origin along this vacuum valley. The wave function can be decomposed as 
dt = fn{f)Xn{u,v\f), with Xn{u,v;r) normalized eigenfunctions of 

H±{u, v; r) = H{u, v, f )\g.=o, see Eq. (|^). It is convenient to use the spherical 
coordinates 

X = f (sin((/3) sin(0), cos((/3) sin(0), cos(0)), (53) 

u = sin^ 0(cos^ 6 + cos^ p sin^ p sin^ 0), v = sin^ 9 cos^ 9 cos p sin p, 

J = J = sin^ 0cos(2(/j) (l — 3sin^ 9 + sin"^ 9[\1 — cos(4(^)]/8) , 

With Xn{9, p\ r) = Xn{u[9, p), v{9, p); r), we define 

3 ^ir/4 i-eicp) 

(xlx') =-^ / dp d9 sin 9 J{9, p)x* {9, p]f)x'{9, p-,r), 

6 J_7r/4 Jo 

OO 

n—O 

where 9(p) implements the constraint X 2 < xs (we can also take 9 S [0,7r] 
and p G [0, 27r], in which case one needs to replace J with | J|/24). We moved 
a factor of 47r to the measure for the f integration, such that /(f) can be 
interpreted as the vacuum-valley wave function (in the S-wave channel). The 
Hamiltonian reduces to 

Xn(-ir-"i4r(^"^r")+<5”'=(En(f)+6f-2))/fc(f), 

n,k,m 

or = S^'^df + Anf), AT^{r) = {Xm\dr\Xn), (55) 


47rf2fif/*(f)/;(f) 


( 54 ) 
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with En{f) the eigenvalues of iJ_L, i.e. H±Xn = En{f)xn- The pairs En{f), 
Xnj can be approximated by the eigenvalues and normalized eigenvectors 
of the truncated reduced Hamiltonian, replacing if (f; L„) with f“^L„(2L„ + 7) 
in Eq. (^o|). From the numerical point of view, for large f the truncation starts 
to become a problem due to the progressive localization of xoj requiring ever 
larger values of L. With L < 20 we can get accurate results for f up to 5. As 
a check, we will also expand f“^i?o(r) and Xo(^! T'-: ici powers of f~^. In the 
adiabatic (or Born-Oppenheimer) approximation which neglects /n> 0 ) Eq. (55) 
becomes 


H'S =Xoi-hr ‘^dff‘^df + Ves{f))fo{f), 

Eeff(f) = ^A^(r) + Eo{f) + 6r“^, (56) 

with A^(r) = — °(^) = (^rXol^fXo)- Unbroken supersymmetry 

requires Ueff(r) = 0. nnP 

An asymptotic analysisocS’El is complicated by the coordinate singulari¬ 
ties. The transverse coordinate is 6 , with exp(—the leading exponential 
factor in xo- Introducing 9 = 9lr^^^, we can expand H± 

^p(o) - - i f 1 

j y 0 l\s J 



82 

o\ 

/o 

0 

o\ 

82 

di 

0 , 

0 

III 

di 

0 

-di 

-82 

0/ 

\ 0 

-di 

0/ 


di = 6 ^dg — 9 %a.n{2tp)d^, 82 = 9 sec{2(j>)dy,, 

A = 9~^dg9^dg + 9~'^ sec{2ip)d^ cos{2ip)dtp. (57) 

There is one exact zero-energy state, xo = ^exp(—i0^), with the proper nor¬ 
malization conveniently expressed as 

Xo = (1 + ^ exp(-iMr^), /i = (2,-2, 2,-4,1, 3)/16. (58) 

2(V7r)'" 

To first order we find {dfXo\dfXo) = 9/(2r^) and Eo{r) = — 33/(2r)^, which 
indeed gives 14ff(f) = 0-|-O(f“®). Note that h ^ but still hN{0)Hf{0)h = 
—2hN{0)h = —2 (see Eqs. ( p^ , ^ ). This can occur because h and differ 
by an element from the (three dimensional) kernel of A^(0). We used H± to 
find Xo, as opposed to N{0)H±, to guarantee (xl^±lxo) vanishes to lowest 
order for any function x{tt,v). The non-trivial kernel makes it particularly 
cumbersome to perform the asymptotic expansion. 
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3.1 Numerical results for Vgff 

Much more convincing, however, are the numerical results depicted in Fig. ^ 
as these “sum” to all orders in f~^. We truncated the matrix {n'\H±\n) to the 
first 420 sphedcal harmonics whose eigenvalues approximate En- Temple’s 
inequality^’N 0 < {H±) — Eq < {[H± — {H±) )/A£', can be used to bound 

the error due to the truncation, with AE = Ei — Eq giving a safe bound. In 
Fig. U this lower bound is indicated by the dots, and the upper bound by the 
drawn lines. The dashed curves (—33/(2r)^ and 9/(2f)^) demonstrate that our 
numerical results perfectly match the asymptotic expansion derived above. 



Figure 4: The energies En(r), for n = 0,1, 2, = ^{SfXoldfXo) and the vacuum valley 

effective potential Fjff (r) (see also the inset). The dashed lines give asymptotic results. The 
dotted curve shows V^fpr) — 6f~^ = i?o(r) + Tho (larger) dots represent the lower 

bound on Eg, using Temple’s inequality. 


We also show Ei{f), to illustrate the gap in the transverse fluctuations. 
The steep dip in this gap around f = 2.5, together with an increasing cou¬ 
pling between the excited states x„, given by {Xn\dfXm), is responsible for 
the breakdown of the adiabatic approximation dramatically illustrated by the 
sudden increase of when approaching r = 2.5. In Fig. ^ we also show E 2 (r) 
to illustrate that the kink in Ei{f) around f = 1.7 is due to an avoided level 
crossing. 

When approaching r = 0, the energies of all the excited states grow as 
due to the non-zero angular momentum, but Eo{f) goes smoothly to zero, as 
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does (9fXo|9fXo)/2- For f ^ 0 one easily finds xo = 2y^j)i = 
and f~^fo(f) will approach a constant (cmp. the behavior of introduced 
in Eq. (^)). Near f = 0 it would be more appropriate to define = 

Eo{f) + as the effective potential, shown in Fig. §as the dotted curve. 

3.2 Groundstate energy and wave function 

The boundary condition to be imposed should be such that /o(fi) goes to 
a constant for large r, and we therefore impose drfoif) = 0 at the bound¬ 
ary of the fundamental domain, f — b = This is equivalent to 

(Xo|5f= 0, where the inner product is at fixed f. Assuming also that 
(Xn|i9f = 0 for all n, we can conclude dr{r^'i') — 0, and the boundary 

condition translates into a condition on the radial momenta, dfji{kb) = 0. 
The b dependence of these momenta and the associated radial matrix elements 
scale with a simple power of b. To obtain the full matrix of the Hamiltonian 
for arbitrary 6, it suffices therefore to calculate momenta and matrix elements 
at 6 = 1. The relevant matrix elements, < > can be computed 

either numerically or following the algorithm developed earlier for the bosonic 
{F = 0) caseEa Needed are the matrix elements with t = 1 {\£' — = 1) and 

t = 4 (|£' — £| = 0,2,4), as well as t = 2, 5 and 8 (with appropriate selec¬ 
tion rules), that occur in determining needed for Temple’s inequality. As 
compared to the bosonic case^ a change in the boundary condition is due to 
the spherical approximation we have used. In reality the boundary is not a 
sphere but a torus, with an appropriately different decomposition of the wave 
function along the vacuum valley. The spherical approximation is justified in 
the supersymmetric case, since fo{r) becomes constant well before we reach 
the boundary. 

In the numerical determination of the groundstate energy for the zero- 
momentum Hamiltonian, there are two reasons this boundary cannot be chosen 
too far from the origin. The first reason is, as seen for Eo{f), that it would 
require too many spherical harmonics to properly localize the wave function in 
the vacuum valley for large b. The second reason is that for increasing b the 
energy gap due to the radial excitations goes to zero, another way of expressing 
the fact that for b —> oo the spectrum becomes continuous down to zero energyH 
To reach 6 = 5 we have used 420 spherical harmonics, and for each up to 20 
radial modes. To keep the size of the matrix manageable, the components 
of the eigenvectors are removed when in absolute value below a threshold 
(typically ~ 10“^), without significantly affecting the accuracy. This process 
of pruning is performed iteratively, increasing the number of radial modes per 
angular state, for those harmonics that stay above the threshold. Together 
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with Temple’s inequality this is extremely efficient to optimize the accuracy 
and achieve numerical control. The number of basis vectors needed increases 
with b, and ranged up to about 3000 to achieve |^ifo| < 0.003 as estimated 
from Temple’s inequality (typically the accuracy of the upper bound is much 
better than this, which we estimate to be l^ifol < 0.0001). 


0 

- 0.02 

- 0.04 


- 0.06 

- 0.08 

- 0.1 


3.5 4 4.5 5 5.5 6 2.5 3 3.5 4 4.5 5 

Figure 5: The groundstate energy Eo(b), using up to 20 radial modes for each of the 420 
harmonics. On the left is shown at an enlarged scale the result obtained with the appropriate 
boundary condition, = 0 (p = 0). The dashed curve is for —8/6^. The same result 

is shown on the right, together with the lower bound from Temple’s inequality (indicated 
by the dots), and in comparison to inappropriate choices of boundary conditions (top curve 
p = 1, lower curve p = —1). 

It is crucial to note that for b finite our boundary condition breaks super¬ 
symmetry. In such a case the groundstate energy need not be positive. To 
illustrate this, and the sensitivity to the boundary conditions, we consider in 
Fig. I the groundstate energy for = 0, with p = —1, 0,1. Indeed, 

p = 0 makes the groundstate energy approach zero most efficiently. It would 
be tempting to conclude this approach is exponential, but our numerical re¬ 
sults rather seem to imply Eo{b) ^ —8/6^, as indicated by the dashed curve in 
this figure. The boundary condition dffo(b) = 0 indeed receives perturbative 
corrections, due to the non-vanishing of dfXoib)- 

However, this is an artifact of the truncation to the zero-momentum modes. 
If we take into account that in the full theory xo al^ involves the non-zero 
momentum modes, the (gauge) symmetry guaranteestj that at the boundary 
of the fundamental domain dfXoib) = 0, and this source of the breaking of su¬ 
persymmetry is absent. The groundstate energy will in this case vanish to all 
orders in perturbation theory. Higher order terms in computing the effective 
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Hamiltonian are required to deal with this, and was the reason for investigat¬ 
ing more closely the determination of the effective Hamiltonian. Corrections 
involving derivatives are manifestations of non-adiabatic behavior, but they 
come from the non-zero momentum modes and can be treated perturbatively 
(even when some of these corrections no longer respect the spherical symme¬ 
try). They were not included here. 



Figure 6: The functions /^(r) (left) and /^(r) — /g(r) (right - note the 40-fold increase in 
scale) extracted from the groundstate wave function ito, satisfying the boundary condition 
9f(f®'I'o)(6) = 0, for b = 4.4, 4,7 and 5.0. We normalized with respect to b = 5, such that 
at r = 2 all agree (indicated by the dot). 


The effect of truncating the Hamiltonian is also clearly seen from the 
behavior of the wave function near the boundary. In Fig. ^ we consider 
the groundstate wave function, satisfying the proper boundary condition 
= 0, and plot p{f) = as well as p{f) - fpf) = 

fniP- With the inner product as defined in Eq. (p^, only involving an¬ 
gular integrations, we recall that fn{f) = o{f)\Xn{r)) (which can be chosen 

to be real), P{f) = f®('I'o(r)|'I'o(f)) and (tpoldto) = 47rr^/^(f)df = 1. In 

the adiabatic, or Born-Oppenheimer approximation f(f) would equal /o('c)- A 
direct measure for the failure of this approximation is given by P{f) — /p (r), 
which as expected deviates from zero when Ei{f) — Eo{f) is small (cmp. Fig. 
but also when we approach the boundary at f = b. In Fig. ^ we show the re¬ 
sults for b = 4.4, 4.7 and 5.0, to illustrate that this deviation decreases with 
increasing b (decreasing coupling, or volume). The results were normalized 
such that for all b, p{f = 2) is equal to its value at 6 = 5. This shows that 
the mismatch between the boundary condition and the truncation of the ef¬ 
fective Hamiltonian does not affect the wave function in the neighborhood of 
the orbifold singularities, where the failure of the adiabatic approximation is 
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non-perturbative. At the same time we see that beyond this neighborhood of 
the orbifold singularities, / will become constant for large 5, compatible with 
a vanishing groundstate energy. 


4 Concluding remarks 

Our analysis shows that, although the orbifold singularities are cause for con¬ 
cern, in the end they do not upset the result for the Witten index. The effective 
Hamiltonian obtained by reducing to the moduli space of flat connections, i.e. 
the vacuum valley parametrized by the abelian zero-momentum modes, requires 
modification due to a singularity in the non-adiabatic behavior at the orbifold 
singularities. One is led to consider the effective Hamiltonian in the full non- 
abelian zero-momentum sector. This removes the singularity. Non-singular 
corrections due to the non-zero momentum modes are still to be included in 
perturbation theory. Supersymmetry is, as usual, expected to keep these per¬ 
turbative corrections in check. In gauge theory the vacuum valley is compact, 
which can be effectively dealt with by imposing boundary conditions in field 
space. Restricting to this fundamental domain is essential, since the non-zero 
momentum modes will give rise to singular non-adiabatic behavior at the other 
orbifold singularities. These are gauge copies of A = 0 and hence outside the 
fundamental domain. The boundary conditions can be argued to preserve the 
supersymmetry if the effective Hamiltonian is constructed to all orders in per¬ 
turbation theory. Despite all similarities, there is an important distinction 
with the problem of the supermembrane, where the truncated Hamiltonian 
is assumed to be exact (apart from approximating SU(oo) by SU(2)) and a 
boundary condition at a finite distance from the origin will always break the 
supersymmetry. This distinction is a subtle, but important one. 

Of course, a numerical analysis can never be entirely conclusive in decid¬ 
ing a theoretical issue that involves the counting of exact zero-energy states. 
Nevertheless, numerical methods do allow us to quantify any non-perturbative 
contributions, which analytically are out of control due to the orbifold sin¬ 
gularities. We have shown that indeed these non-perturbative effects do not 
contribute to the vacuum energy, which thus remains zero. In this paper 
we have only considered SU(2), to illustrate how to go beyond the adiabatic 
approximation. There is no fundamental obstacle to consider other groups, 
including dealing with the newly found disconnected vacuum components, but 
technically this will be much more demanding. 

An additional motivation to push ahead with this approach was that our 
methods and results may also be relevant for more general situations in which 
the zero-momentum Hamiltonian (in its truncated form) seems to have a role 
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to play. To this purpose we carefully documepled our computer code, and 
make it available through the World Wide WebO 
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Appendix 

We can also diagonalize Hf{c) by solving the system S — EfS and V = EfV 
(see Eq. (^). Assuming Ef ^ 0, we solve EfV = EjV hy replacing EfS, as it 
appears in EfV, hy S (which does not contain S). This gives a linear system 
of equations for V, which is however quadratic in Ef, 


MtVt = 2ctc\vt + crclVt + cic\Vf + EfS^.^e^ 


abd 

'3 


c-M = EjVi 


(59) 


For c“ = diag(a:i, X 2 , xs) (no summations over repeated indices) 

— 2iXiXk5jia^th J- Xikf^SdiSkb T t^k^ki^ab E Ef 6}zij£abj^j , (1^^) 


which splits in the 3x3 block AI§® = Al|- = 2xl5kiE2xiXk+Ef{J\- Xj)/{xiXk) 
and three 2x2 blocks (forming a triplet under the Xi permutation symmetry) 
Mf = M^l = XiXk + EfXj - EfXjSik , (a, b are fixed by requiring SkajSibj 7 ^ 0, 
which also constrains i,k ^ j). Only AIq will give rise to invariant states, 


det(AIo - Ej) = -{E^f - Ac^Ef - 16detc)(E;^ - c^Ef + 2detc) = 0, (61) 

agreeing with Eq. (p^ ) . Also note that det (AI^ — Ej) = Ef{Ej — c^Ef + 2 det c) 
for all j, and its resulting energies {Ef ^0) agree with the energie^f the non¬ 
invariant two-gluino states constructed from the one-gluino statesEa 
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